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Abstract 

We transform the two-matrix model, studied by P.Di Francesco and al. in 
into a normal one-matrix model by identifying a "formal" integral used by 
these authors as a proper integral. We show also, using their method, that the 
results obtained f or the resolvent and the density are not reliable. 



1 THE MODEL 

In a recent paper P. Di Francesco, B. Eynard and E. Guitter |[TJ discuss a model of 
two n x n hermitean matrices M and R with a partition function 

Z{p,q,g;N) = J dMdRexp{-NTr[p\og(l - M) + q log(l - R) + gMR]} (1) 

where g is later replaced by 

» = i (2 > 

and the strong coupling limit is of primary interest: t — > 0. Applying the Itzykson- 
Zuber integral identity, the partition function is transformed into an integral over the 
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eigenvalues {mj}j=i v . n and {ri}i = i t .. n of M and R respectively. 

/n 
A(m)A(r) n[exp{-iV[plog(l - m*) (3) 
i=i 

+glog(l - n) + gmiTiiydridmi] 

where A is the Vandermonde determinant. 
However, an integral such as 

/ dxdyexp{-N[p\og(l - x) + <?log(l - y) + 0a?y]}j? n (a;)p TO (y) (4) 

where p n {p m ) are polynomials of degree n(m) does not exist. 

In later parts of the paper the authors want to ascribe a meaning to such integrals 
by the "formal integral" (see eq. (4.2)) 

[ dxdy x a y? exp{-^-xy] = a\ S a(3 (^-) a+1 (5) 
Jsr 2 t TV 

Special cases of such "formal integral" are already used earlier in the text. It is, 
however, easy to see that the "formal integral" is a proper integral 

2 [ dxdy w a w f3 exp[-—ww] = a\ 5 a g (4)" +1 ( 6 ) 

JE 2 t N 

(w = x + iy , w = x — iy) 

It is also easy to go back and replace all "formal integrals" by proper integrals, then 
for (|3]) we obtain 

„ n 

Z(p,q,g;N)~ j \ A(w) \ 2 n[exp{-JV[plog(l - w t ) (7) 

i— 1 

+glog(l - Wi) + gWiWi]}dxidyi] 
where we have to restrict the Wi to, say 

<| Wi |< 1 (8) 

and 

Re{Np) < l,Re(Nq) < 1 (9) 

if we want to give an analytic meaning to Z in the variable t = g^ 1 , but we can 
integrate over the whole complex plane if we are interested only in the formal power 
expansion in t for t — > 0. Finally, without referring to the Itzykson-Zuber formula we 
can integrate over U (n) to get 

Z(p,q,g;N) = T J dMdM^ exp{-A^Tr[plog(l-M) + glog(l-M t ) + c/MM t ]} (10) 

where M is a normal matrix 

[M,M*] = Q (11) 

Thus the model actually evaluated is not a hermitean two-matrix model but a normal 
one-matrix model. 
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2 THE SADDLE POINT EQUATION 

The partition function is calculated in the limit iV — > oo so that n = zN and z is 
kept fixed. In this limit the partition function is approximated by (see eq.(4.20) of 

0) 

/•CO 

Z = N n dai...da n exp{N 2 S(ai..a n ,p,q, z)\ (12) 
Jo 



with 



1 n 

S = —J2 I (ai + P - z)(\og(ai + p - z) - 1) 



N i=1 

+ (aii + q- z)(log(a!j + q - z) - 1) - 2a i (log(Q; i ) - 1) 
1 N t 

+ _iog( 7 (Q^v+i,-))+Q*iog(-: 

+ ^E lo g(l -«i I) 

- [ ds[(p - s)(log(p - s) - 1) + (q - a)(log(g - s) - 1) 

JO 

+s(log(a)-l)+slog(t)] (13) 

where 7(0;, £) is the incomplete 7-function p. In the limit N 00 the term contain- 
ing the incomplete 7-function is calculated as 

where £ = at and so that the derivative of the r.h.s is continuous at £ = 1. The 
saddle point equation is then 

2P(dfi4®- = Q(l-at)\og(t {a + p - Z)ia + q - Z) ) 
J p — a v a ' 



f3 — a 

(a + p — z) (a + q — z) • 



+ e(at-l)]os( ^ P ~ Z ^ q - Z) ) (15) 
valid for a in the support of the density p and with the normalisation 

d/3p(P) = z (16) 



We mention that the free energy also has an additional term compared with the 
expression obtained in the equation (4.29) of Q 

/■OO 1 1 /"OO 

td t f = / dppptf) - ^ + 7 / d/3p(p) (17) 
jo It Jt- 1 
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We solve the saddle point equation by making the ansatz 

p(a) 
Pi{a) > 
z\ + z 2 

and 



supp(pi) 
supp{p 2 ) 



Pi (a) + Ma) 

d(3pi((3) = Zi 



< 7b 72 > 

< 73,74 > 



(18) 
(19) 
(20) 



(21) 
(22) 



In order to solve (|T5[) we have to make sure that 



72 < - < 73 



Next we introduce the resolvent 



u>i(a) = I df3 



. PiiP) 
a — (3 



so that ( |i5|) induces 

u>i(a + iQ) + u)\(a — iO) 



log (t 



uj 2 {a 



+ z0) +u 2 (a-i0) = -log( 



(a + p - z) (a + q - z) 
a 

(a + p - z)(a + q - z) 
a 2 



and 



lim auAa) 

a— >oo 



d(3pi({3) = Zi 



(23) 



(24) 



(25) 
(26) 

(27) 



We shall use the method of 0, and show that their result for u>i, p\ is not reliable, in 
spite of the fact that the saddle point equation (|25|) is the same as the one obtained 
in [0]. For simplicity we set Z\ = z, z 2 = from now on. 

We introduce two parameters r, 5 and three hyperbolic angles (assumed to be all 
> 0)$i,$ 2 ,$3 by 



a = z — r — 25cosh($) 

p = r + 2<5cosh( ( I ) i) 

q = r + 28 cosh($ 2 ) 

z = r + 25 cosh($ 3 ) 



Then 



a 



-(a + p — z)(a + q — z) = 5t- 



4(T 1 2 -T 2 )(T|-T 2 )(l-r|) 
l_ T 2 )(1 _ r 2 )(1 _ T2)(T 2_ T2) 



(28) 
(29) 
(30) 
(31) 



(32) 
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where 

T = tanh(-) (33) 
Ti = tanh(^)>0 (34) 

with 

71 = z — r — 26 , 72 = z — r + 26 

Along the cut < 71,72 > the parameter $ is assumed to vary as +iip (0 < cp < tt) 
for Im(a) \ and as — iip for Im(a) /* 0. On the real axis outside the cut, $ is 
necessarily so that 

$ 

T = tanh(— ) — > —1 for a — > ±00 (35) 

Taking into account Q3"3) ) and ( |33| ) only one ansatz for Wi(a) is possible which has 
only the cut at < 71, 72 >, namely 

1 1 1 1 2(r 1 -r)(r 2 -r)(i + r 3 ) ^ 



provided we take 



(i-71)(i-r 2 )(i + r 3 ) 

(1 + 71)(1 + T 2 )(1-T 3 ) 
= u 1 u 2 u^ 1 (37) 



where 



l+Ti 

1-T 

u = — = e 

1 +T 



Finally we obtain for the density 

Pl (a) = ^(<P- 2(*i(p) + *2(y) - * 3 (*>))) (3? 

where 



/ x ,tan(^ N 
ty(a) = arctan(- 



We can show that 



ckp26san(<p)px(J3(<p)) 

6[u 3 -u x - u 2 ] 
lim aui(a) (39) 
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as desired. 

Positivity of the density function (^) is achieved if and only if 

T x + T 2 - T 3 - 1 > (40) 

This constraint follows from the observation that the neighborhood of 7r in the variable 
tp is critical, i.e. negative values of p\ appear here first. Namely, if we write 

\ = \-X (41) 



we have then 



and consequently 
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